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1 Introduction 



The point of view developed in this article is, given some continuous differential 
equation, to find its discretizations which preserve some global property, namely 
the explicit linearizability or more generally the discrete Painleve property, a notion 
defined section ^| together with its group of invariance. From this point of view, 
we will never have to address the question of finding a continuous limit, and conse- 
quently we will always write the discrete equations in a form as close as possible to 
the canonical form of the continuous limit, which is well established [0J |2|. 

Discrete equations can be considered as functional equations linking the values 
taken by some field variable u at a finite number iV+1 of points, either arithmetically 
consecutive : x + kh, k — ko = 0, 1, . . . , N, or geometrically consecutive : xq k , k — k = 
0, 1, . . . , N, where h or q is the lattice stepsize, assumed to lay in some neighborhood 
of, respectively, or 1, and ko is just some convenient origin. The integer N is called 
the order of the equation. They have received a decisive impulse from statistical 
physics in 1990 with the discovery of a discrete analog of the first Painleve equation 

(pi) i a § 

E = -(u-2u + u)/h 2 + 2(u + u + u)u + x = (1) 

and of the second Painleve equation (P2) || [7j (in the particular case a = 0) 

E = -{u-2u + u)/h 2 + (u + u)u 2 + xu + a = 0. (2) 

The above short notation 

u = u(x), u = u(x + h), u = u(x — h), u = u(x + 2h), ... (3) 

is adopted throughout the article. These two equations clearly admit as continuous 
limits (h -> 0) (PI) and (P2) 

E = -d 2 u/dx 2 + 6u 2 + x = 0, (4) 
E = -d 2 u/dx 2 + 2u 3 + xu + a = 0. (5) 



In fact the equation (|I|) had been written down by many authors j8], |9], [L0], [TT|, [13 
but ref. 0, |, |] were the first to notice the continuous limit to (PI). 

The reason why these two discrete equations deserve the name of discrete Painleve 
equations, in short d-(Pl) and d-(P2), is that they possess the discrete Painleve 
property. Indeed, both admit a discrete Lax pair (A,B,z,ip,h) JOJ 

i[)(x + h/2) = A%l)(x - h/2), d z ij(x - h/2) = B^{x - h/2). (6) 

The matricial Lax pair of the d-(P2) (H) is P], BO] 



fl/z hu\ z — z 1 , „ /7 9n 

A ={L z)> t = ^2h-> Z = l + ht + ° {h) 

hzB = (-2uu + At 2 - {x - h/2)) (I °A + (z + l/z)(u-u)/h^ l J 

-(2t(u + u) + aW? I) (7) 



d z A + AB — BA = 2z ( °, \ I E, 



1 
1 

-1 
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while the matricial form of that of the d-(Pl) ([!]) is 



A = ( 1 ~f Z h f U _~ h2 f ),t = z- h 2 z*/2, z = t + (t 2 /2)h 2 + 0{h*) 
2(1 - h 2 z){u - u)/h -4uu - 2{x - h/2) + At(u +u)- 16t 2 ' 



B 



2(u + u) + 8t -2(1 - h 2 z){u -u)/h 



h-\d z A + AB-BA) = 2Q ) 

a new result to our knowledge. 

Just like its continuous counterpart, the discrete Painleve test is the set of all 
methods one can imagine to build necessary conditions for a given discrete equation 
to possess the discrete Painleve property. At the present time, only one such method 
is known, namely the singularity confinement method of Grammaticos et al. ||15||. We 



propose here a new method, based on the perturbation theorem of Poincare, which 
is applicable to differential systems in an arbitrary number of independent variables, 
whether discrete, continuous or even mixed discrete-continuous. As an example, we 
apply it to a qualitative discretization of (PI) and isolate three candidates d-(Pl). 

The paper is organized as follows. The discrete Painleve property (PP) and 
its group of invariance are defined in section |2], and a first set of basic rules of 
discretization is given. Section |3| explains the new method for the discrete Painleve 
test, introduces a new notion, the depth of a discrete equation, and applies the 
method to qualitative candidates d-(Pl) and d-(P2). 



The discrete Painleve property, its group of in- 
variance, basic rules of discretization 



The (continuous) Painleve property is defined [I(J as the absence of movable 
critical singularities in the general solution of a differential equation 

Va; : E(x, u, du/dx, . . . , d N u/dx N ) = 0, (9) 

where a singularity is said movable (as opposed to fixed) if its location in the complex 
plane of x depends on the initial conditions, and critical if the solution is multivalued 
around it. For shortness, following Bureau [pj]], we will use the terms "stability" for 
PP, "stable" or "unstable" for an equation with or without the PP. 
The PP is invariant under the group of birational transformations 

(u, x) -+(U,X): u = r(x, U, dU/dX, . . . , d N ~ 1 U/dX N - 1 ) = 0, x = H(X)(,10) 
(U, X) -+(u,x): U = R(X, u, du/dx, . . . , d^u/dx"' 1 ) = 0, X = f (x),(ll) 

(r and R rational in U, u and their derivatives, analytic in x,X). An easier to 
manage subgroup is made of the homographic transformations 

(u,x)^(U,X): u = ^-^ X = ax), ad-bc^0 (12) 

cU + d 

where (a, b, c,d,£) are arbitrary analytic functions of x. In his classification of second 
order first degree equations, Gambier || has found respectively twenty-four and fifty 
equivalence classes for these two groups (with minor later corrections [I8[ |19[). 
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In the discrete case, let us consider equations 

Vx V/i : E(x, h, {u(x + kh), k — k = 0, . . . , N}) = (13) 
VxVg : E(x,q,{u(xq k ), k — k = 0, . . . , N}) = (14) 

algebraic in the values of the field variable, with coefficients analytic in x and the 
stepsize. It should be noted that u is a function of two variables, x and the stepsize. 
A natural definition for the discrete Painleve property, which seems to have never 
been given before, is the following. 

Definition. A discrete equation is said to possess the discrete Painleve property 
if and only if there exists a neighborhood of h = (resp. q — 1) at every point of 
which the general solution x —>■ u(x,h) (resp. x — > u(x,q)) has no movable critical 
singularities. 

Remarks. 

1. The definition reduces to that of the continuous PP in the continuous limit. 

2. This definition immediately extends to equations in an arbitrary number of 
independent variables, discrete or continuous, the extension starting then from 
the definition of the PP suited to partial differential equations (PDEs), which 
we do not remind here since this is not our subject. 



The discrete PP is invariant under the discrete analog of (|Tl"D, which is the group 
of nonlocal discrete birational transformations 

u = r(x, h or q, U, U, U, . . .), 

U = R(X, H or Q, u,u, u, . . .), X = £(x, h or q), H = r](h), Q = n(q), (15) 

(r and R rational in U,U,U_, . . . ,u,u,u, . . ., analytic in x and the stepsize, £, r], k 
analytic). There exist two discrete analogs of the subgroup (|i~2"D, and both may 
be useful to establish the discrete equivalent of the classification of Gambier. The 
first one is the group of nonlocal transformations (|T5|) which in the continuous limit 
reduce to the homographic transformations (|T2"D, where (a,b,c,d,l;) are arbitrary 
analytic functions of x and of the stepsize. The second one is the group of local 
homographic transformations (r and R homographic in U and u, independent of 
U,U, . . . ,u,u, . . ., analytic in x and the stepsize, £, rj, k analytic). 

Remark. The first subgroup seems more useful, although it does not contain the 
transformation u = h h U, k G Z. 

Let us now give some basic rules for discretizing a given continuous equation (|9]) 
into either a difference equation (|T3|) or a g-difference equation (|TJ). 

The question of discretization is well known in numerical analysis, where one 
looks for a scheme of discretization which maximizes the order, called scheme order, 



of the remainder of the expansion of the left-hand side of (|Tj) in a Taylor series of 
h around the center of the N + 1 points. A scheme of discretization is said exact iff 
it has an infinite order, like the one for the linear ODE u {N) = 

N 

J2(-l) j C' N u(x + jh)=0. (16) 

3=0 

Contrary to numerical analysis, only interested in a local integration, we require 
the scheme of discretization to preserve the differential order N, an essential element 
for a global knowledge of the solution : every discretization must involve iV+1 points. 
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Just like in the continuous case, there exists another important element concern- 
ing discretization rules. 

Definition. The degree of a discrete equation is the highest of the two polynomial 
degrees of the LHS E of the equation in u(x) and u(x + Nh), or u(x) and u(xq N ), 
where E is assumed polynomial in the N + 1 variables u(x + kh) or u(xq k ). 

Let us have a closer look at the degree before and after a discretization. 

Consider first a second order, first degree equation with a single valued general 
solution 

u = ( Cl x + c 2 ) 2 , uu" - (1/2V 2 = 0, (17) 

and its exact discretization 

u-2u + u , . (u-u)(u-u) 2 fu — 2u + u\ 2 n , 10 . 
u (1/2P ± >--{l/%)h 2 [^ J =0. (18) 

Although the exact discretization does not conserve the degree, there exists a second 
order scheme which conserves it : this is the one defined by dropping the third term 
in© 

u U-2u + u _ m (u-uXu-u) ^ Q (ig) 

What is remarkable is that both schemes are explicitly linearizable into ip — 2if) + ip = 
by the respective transformations 

u = ^\ _ (20) 

u= (l/4)[^(V + ^) + V 2 + ^]- (21) 

Consider next a first order, first degree ODE with a multivalued general solution 
u = ( x - x )~ 1/2 , u' + (l/2)w 3 = 0, (22) 
and its exact discretization 

u — u u 2 u 2 



, + = 0. (23) 

h u+u v ' 

There exists no discretization scheme conserving the degree one of the continuous 
equation, and the reason is here the multivaluedness. 

We therefore conjecture : given an algebraic differential equation with the PP, 
there exists a discretization scheme of order two which conserves the degree. 

Hence the additional rule, restricted to algebraic differential equations with the 
PP : every second order scheme must also conserve the degree. 



3 The discrete Painleve test 

In the continuous case, all the methods of the Painleve test, without exception, 
are based on two theorems and only two, namely the existence theorem of Cauchy 
and the theorem of perturbations of Poincare |21[ . This is explained in detail in the 
lecture notes of a Chamonix school p2 |. 
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The discrete case is an arithmetics problem, hence very difficult. But a usual 
approach in arithmetics is to convert the problem into one of analysis, for which all 
the tools dealing with the notion of limit are available. This is appropriate here, in 
order to make applicable the perturbation theorem of Poincare. 

Consider an arbitrary discrete equation (0), also depending on some parame- 
ters a, and let (x, h, u, a) — > (X, H, U, A, e) be an arbitrary perturbation admissible 
by the theorem of Poincare (which excludes any nonanalyticity, like e 1//5 , for the 
perturbed variables). A necessary condition is that the limit e = possesses the PP 
(discrete or continuous, this does not matter). 

In order to illustrate the method, let us discretize the two equations (PI) and 
(P2) by a second order scheme, using the rules previously stated. We will detail the 
case of (PI) and just give the result for (P2). 

The candidates d-(Pl) and d-(P2) are as follows (we choose them rather sim- 
plified to alleviate the presentation) 

E = -(u-2u + u)h~ 2 + 3Ai(m + m)m + 6A 2 m 2 + 6\ 3 uu + g = 0. (24) 
E = -(u - 2u + u)h~ 2 + Ai(m + u)u 2 + 2\ 2 u 3 + 2\ 3 uuu 

+x(fj, 1 (u + u)/2 + /j, 2 u) + a = 0, (25) 

with Yl Afc = 1, Yl fJ'k — 1) a a constant and g an unspecified function of x. 

The test will generate necessary conditions on {X^^g). One must find at least 
the following solutions, where the equations have a Lax pair : for d-(Pl), g = x 
and X = (2/3,1/3,0); for d-(P2), X = (1,0,0), jt = (0,1). One may also find 
the following solutions, isolated by the singularity confinement method although 
no Lax pair is yet known : for d-(Pl), g = x and A = (1,0,0) (ref. ||15||), A = 
(1/2,1/4,1/4) (ref.H eq. (5.5)). 



The method which we propose is a perturbation of the continuous limit entirely 
analogous to the Fuchsian [^3| or nonFuchsian |24| perturbative method of the con- 



tinuous case. This perturbation is defined by an expansion of u as a Taylor series 
in e, where e is nothing else than the lattice stepsize 



oo 

x unchanged, h = e, q = e £ , u = ^2 £ n u^ n \ a = analytic (A,e). (26) 



n=0 



It generates an infinite sequence of (continuous) differential equations E^ = 
defined by 



E = J2 e n E {n) (27) 

n=0 



whose first one n = is the "continuous limit". The next ones n > 1, which are 
linear inhomogeneous, have the same homogeneous part E^ u^ n ' = independent 
of n, defined by the derivative of the equation of the continuous limit, while their 
inhomogeneous part ( "right-hand side" ) comes at the same time from the non- 
linearities and the discretization. 

The only tiny difference with the continuous perturbative method lies in the 
explicit dependence on e of the unperturbed equation. This does not affect the 
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rest of the method and, depending on the nature, Fuchsian or nonFuchsian, of the 
linearized equation E^ l > = at a singulier point of one then applies, without 
any other change, either the Fuchsian perturbative method of the continuous case 
(Ref. [|23| p. 42) or the nonFuchsian perturbative method of the continuous case 
(Ref. [TJ|] p. 342). We adopt the notation used in these two articles. 

As a tutorial introduction, let us first handle the Euler scheme for the Bernoulli 
equation 

E = {u-u)/h + u 2 = (28) 

i.e. the logistic map of Verhulst, a paradigm of chaotic behaviour. Let us expand 
the terms of ( |28| ) according to the perturbation ( |26| ) up to an order in e sufficient 



to build the first equation = beyond the continuous limit E^ = 



u = «(°) + u w e + 0(e 2 ) (29) 
u 2 = u^ 2 +2u^u^e + 0(e 2 ) (30) 
u = u + u'h + {l/2)u"h 2 + 0{h 3 ) (31) 

?Lz?L = tt («)' + ( M W' + (i/2)/)") e + 0( £ 2 ). (32) 

The equations of orders n = et n = 1 are written as 

E<® = u W+u^ 2 =0 (33) 
E (i) = E (oy u (i) + (1/2 ) u (o)" = o, £(<>)' =dx + 2w (0). (34) 

Their general solution is 

u = x~ 1 ^X = x ~ x o, xq arbitrary (35) 
= w_iX~ 2 — x " 2 Log-?/', ip = x — xq, ul\ arbitrary, (36) 

and the movable logarithm proves the instability as soon as order n = 1, at the 
Fuchs index % — — 1. 

Definition. One calls depth of a discrete equation the rank n of the first Taylor 
coefficient of u in series of the stepsize not to possess the continuous PP. 

The highest the depth, the largest the choice of the stepsize for a numerical 
integration. Equal to one only for the Verhulst map, the depth is infinite for any 
stable equation. 

Remarks. 



1. The depth is bounded from above by the smallest of the perturbative orders n 
at which a movable logarithm appears when one runs over the several families 
of movable singularities. 

2. The only restriction on is not to be what is called a singular solution (not 
obtainable from the general solution by assigning values to the arbitrary data), 
i.e. it can be either the general solution (as above) or a particular one, it can 
also be either global (as above) or local (Laurent series). 

3. The derivative of E^ at having iV as differential order since the singular 
solutions are excluded, the sum + u^e is a local representation of the 
general solution in the neighborhood of (x, £ ) — (0,0). One must of course at 
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the order n — 1 include in vy) all the arbitrary data missing in m' ', and only 
them. Here, no such arbitrary need be included, and one can set to zero the 
term u^\x~ 2 , i.e. the partial derivative of with respect to the arbitrary 
data Xq. 



For the example (23), one obtains 

E (o) = _ u (o)" + 6u (of +g ( x) = Q (37) 

EW = -dl + Uu®, (38) 
R {1) = (39) 
R (2) = 6M (i) 2 _ 6 A3 M ( )%3(A 1 + 2A3)A (0 )''-(l/12) M (°)'''', (40) 

etc. At order n = (the "continuous limit"), Painleve has found the necessary 
condition g" = 0. Let us choose g(x) = x + b and for the local general solution 
of fl37|) 

u<°> = x- 2 (l-(&/10)x 4 -(l/6)x 5 + 4 0) X 6 + --.) (41) 

in which a;o and Mq '' are arbitrary. This is then the Fuchsian perturbative method 
which must be used 

n > 1 : W W = £ uf X j+P , = 0, 4" } = 0, = £ ^] n) X J+9 , (42) 



in order to generate the necessary conditions 

n > 1 : Q, W = = 0, i G {-1,6}. (43) 

The first nonzero conditions are 

= 3(-A 2 + 3A^ + 10A 3 + 21A 2 A 3 - 60A§)/50 = (44) 
Q {6 \ = 9(-A 2 + 2X 2 2 + 3X 3 2 - 19A 3 + 28A 2 A 3 

+11A^A 3 + 58A^ - 7X 2 \l - 39Alj)/5 = (45) 

Q< 8) = 6/4 = (46) 

gPJ) = 6/5 = (47) 

Q^ 10) - 4°V 5 = (48) 

= 4 0) / 6 = (49) 

Q$*> = b 2 U = (50) 

in which the fk's, all different, represent polynomials of degree k in (Ai, A 2 , A3). 
The polynomials in the variables (Ai, A 2 , A 3 ) are generated by the ideal 

Ai + A 2 + A 3 - 1, A 3 (A 2 - A 3 ), A 3 (Ai - 2A 2 ), A 2 (A X - 2A 2 ), (51) 

and this ideal possesses the only three zeroes : 

(Ai,A 2 ,A 3 ) = (2/3,1/3,0), (1,0,0), (1/2,1/4,1/4). (52) 
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This result is obtained at (n, i) = (6, —1), but one must proceed to (n, i) = (8, 6) to 
get rid of three irrational extraneous solutions; we have checked the vanishing of all 
no- log conditions up to (n, i) = (22,6). 

The first value A = (2/3,1/3,0) (case a = 1 in Ref. ]n|) corresponds to the 
d-(Pl) (HD with a Lax pair found in quantum gravity, the condition is then sufficient. 

The second value (1,0,0) (case a = in Ref. ]n|) corresponds to a candidate 
d-(Pl) which admits the second order matricial Lax pair 

A _(a x 2h(f 2 u - / 3 
A ~\Q -Oi 
# 2jl = 0, B 1A = -B 2j2 = 4J(h~ 3 ai /f 2 + #/(2J/ 2 )), 
S 1>2 = J(-12uu - 2(x - /i/2) + 4/T 2 (u + u) - 8h~ 2 f 3 /f 2 ) 
+ half's- {hi f 2 )f 2 ) 

d z A + AB-BA=(° Q " 4 Q lJ )i?, (53) 

in which J, / 2 , f 3 are arbitrary functions of z, and ai an arbitrary constant; however, 
we have not succeeded in finding a continuous limit to this Lax pair. 

The third value (1/2, 1/4, 1/4) defines an equation E(u) = equivalent to that 
E(v) = for (1, 0, 0) under the discrete birational transformation ^6] 



2v(x) = u(x + h/2) + u(x - h/2) 

2u(x) = -3h 2 v(x - h/2)v(x + h/2) + v(x - h/2) + v(x + h/2) - /i 2 x/2.(54) 



As to the scheme ( p5|) for d-(P2), it admits the unique solution A = (1,0,0), 
Jt = (0,1), obtained, for the family u' ' ~ ±X _1 with Fuchs indices (—1,4), at 
(n,i) = (4,4). This is equ. (0), so the condition is also sufficient. 

Remarks. 

1. These four solutions are simple rational numbers, and moreover they are con- 
vex : Wk : < A fc < 1, < < 1. 



Following Painleve p7|| , one should in fact search for the "complete equation", 
i.e. for all the admissible nondominant terms which can be added to the can- 
didate d-(Pl) ( |24"|) without destroying the PP. This will be done in future 
work. 



3.1 Discussion 



The example fl2"l|) has been chosen for its simplicity but it already allows for 
some comparison between the singularity confinement method and the perturbation 
of the continuous limit. 

When applied to (24]), the singularity confinement method finds the same nec- 
essary conditions ((52]) , as we have checked (this was checked earlier with different 
types of equations |15|, |26f). So the most important point is, again for this particular 
example, the identity of the set of necessary conditions generated by both methods. 
Other examples are currently under investigation |32 , |28 | to detect situations where 
the two methods would produce complementary, not identical results. 
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We have no explanation for the apparent identity of the results produced by the 
two methods, but a reason could be the existence of a perturbation of Poincare iden- 
tical to the method of confinement, and whose no-log conditions would be identical 
to the confinement conditions. Despite a thorough search for such a perturbation, 
we have not found one. 

Since it only relies on the existence of a continuous limit, our method can be 
extended without difficulty to equations in an arbitrary number of dependent or 
independent variables, whether the equations be discrete or mixed continuous and 
discrete. 

The singularity confinement method has also been extended to such situations 



2!|, |30[; however, in the case of m discrete independent variables, one must check in 
addition that the result of the iteration is independent of the path followed on the 
m— th dimensional lattice. 

As pointed out by a referee, the present method is also applicable to a discrete 
equation admitting several continuous limits. In such a case, the test must evidently 
be performed around every continuous limit. 



4 Conclusion 

We have tried to give some precise definitions and guidelines to make more 
systematic the search for discrete analogs of differential equations integrable in the 
sense of Painleve. 

The discrete Painleve test is now made of two main methods : the singularity 
confinement method, which seems intrinsically discrete, and the present one, which 
assumes an underlying continuous limit. Both points of view (discrete first, contin- 
uous first) are useful since they generate identical results on a sample of equations. 

Rules of discretization for Lax pairs will, in future work, be established and 
systematically applied to the numerous d-(Pn) or q— (Pn) candidates which still 
lack a discrete Lax pair. 

Finally, let us mention that the celebrated equation of class III of Chazy |31] 



v!" — 2uu" + 3m' 2 = 0, has recently been discretized [^2j following the methods of 



this paper : the no-log conditions are entirely obtained at order n = 2, and then 
they are all satisfied at least up to n = 16. However, this candidate d-(Chazy) has 
not yet been integrated. 
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